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We use the finite-size density-matrix-renormalization-group (FSDMRG) method to obtain the 
phase diagram of the one-dimensional (d = 1) extended Bose-Hubbard model for density p = 1 in 
the U — V plane, where U and V are, respectively, onsite and nearest-neighbor interactions. The 
phase diagram comprises three phases: Superfluid (SF), Mott Insulator (MI) and Mass Density 
Wave (MDW). For small values of U and V, we get a reentrant SF-MI-SF phase transition. For 
intermediate values of interactions the SF phase is sandwiched between MI and MDW phases with 
continuous SF-MI and SF-MDW transitions. We show, by a detailed finite-size scaling analysis, that 
the MI-SF transition is of Kosterlitz-Thouless (KT) type whereas the MDW-SF transition has both 
KT and two-dimensional-Ising characters. For large values of U and V we get a direct, first-order, 
MI-MDW transition. The MI-SF, MDW-SF and MI-MDW phase boundaries join at a bicritical 
point at (17, V) = (8.5 ± 0.05, 4.75 ± 0.05). 

PACS numbers: 05.30Jp,67.40Db,73.43Nq 
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I. INTRODUCTION 



model whose Hamiltonian is 
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+V (!) 

<i,j> 

The first term in Eq. Q represents the kinetic energy- 
associated with the hopping of bosons from site i to 
its nearest-neighbor site j with amplitude t; a\ (dj) is 
the boson creation (annihilation) operator at site i and 
hi = a\di is the associated number operator; onsite U 
and nearest-neighbor V interactions are represented, re- 
spectively, by the second and third terms and are pos- 
itive since they are repulsive. We restrict ourselves to 
the physically relevant region V < U and set the scale of 
energies by using t = 1. 

Model (JTJ has been studied by a number of authors^ 
in the case V = 0, i.e., in the absence of nearest- 
neighbor interactions; at zero temperature (T = 0) it 
has been shown to have a superfluid phase if the mean 
number of bosons per site p is not an integer; how- 
ever, for integer densities, it shows a superfluid (SF) to 
Mott insulator (MI) transition. This SF-MI transition 
is of the Kosterlitz-Thouless (KT) typeSi in one dimen- 

: i c n 2.7.8.9.10.11 

In the limit U — * oo model maps onto the spin- 
i XXZ model if the mean number of bosons per site 
p = 1/2. Every site can now have only two possible 
states, namely, a state with no boson and another with 
one boson. We represent these as | 0) and | 1), re- 
spectively, and make the identification | 0) — > |J.) and 

1) — *|T), where \\) and ||) are, respectively, spin-i 
down and up states. Now, by using the transformations 
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The study of quantum phase transitions in systems 
of interacting bosons is an exciting area with a fruitful 
interplay between theoryi2i2i2i£i§iL2iiLiiI, numerical sim- 
ulationsAi*i2iiii4, and experiments. A variety of experi- 
mental systems have been studied: liquid 4 He in porous 
media like vycor or aerogel^; microfabricated Josephson- 
junction array s 16 ' 17 ; the disorder-driven superconductor- 
insulator transition in thin films of superconducting ma- 
terials like bismuthA^; flux lines in type-II superconduc- 
tors pinned by columnar defects aligned with an exter- 
nal magnetic fieldi^; and best from the point of view of 
comparing theory with experiments, atoms trapped in 
optical-lattice potentials. In a system where the num- 
ber of atoms per site is an integer, Greiner et aim* have 
observed a superfluid-Mott insulator transition for 87 Rb 
atoms, trapped in a three-dimensional optical-lattice po- 
tential, by changing the strength of the onsite potential. 
Experiments in such optical lattices have several advan- 
tages over their condensed-matter counterparts, includ- 
ing precise knowledge of the underlying microscopic mod- 
els2i, the possibility of controlling parameters in the ef- 
fective lattice Hamiltonians, and the absence of disorder. 
The recent probable observation of a supersolid helium 
phase 22 has given a further fillip to this area. Even in 
the absence of disorder these systems can show a variety 
of phases like Superfluid (SF), Mott Insulator (MI), and 
Mass Density Wave (MDW) [or Charge Density Wave 
(CDW) if the bosons are charged]. The simplest model 
that can show these phases is the extended Bose-Hubbard 
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at = Sf, a,i = S i , and n, = (Sf - |), the model (JTJ 
maps onto the spin-^ XXZ model with the Hamiltonian 

Hxxz = -2t ( S f S * + S ! S j) + V 2 S ' S '> (2) 

<i,j> <i,3> 

where we have suppressed constant terms. This model 
has been solved exactlji^i and shows a KT-type transi- 
tion from XY to Ising ordering at V = 2t. The bosonic 
analogs of XY and Ising phases are, respectively, SF and 
MDW phases. 

If p = 1 and t = 0, it is easy to see that model 
has a first-order, MI-MDW transition at U = 2V. Large 
values of U favor the MI phase whereas large values of V 
favor the MDW phase. 

Recently Kiihner et aiS studied model (JTJ in one di- 
mension by the using a finite-size, density-matrix renor- 
malization group 2 ^ (FSDMRG) and showed that, for 
V = 0.4, it has a continuous SF-MDW transition for den- 
sity p — \ and a continuous SF-MI transition for p = 1. 
Niyaz et a& have studied this model in one dimension by 
a Monte-Carlo method and obtained its phase diagram 
in the (U, V) plane for p = 1. They obtain SF, MI and 
MDW phases in model (JTJ and continuous SF-MI, SF- 
MDW, and MI-MDW transitions but conjecture that, at 
large U, the MI-MDW transition should be first order. 
The study of Ref. [2|| has obtained a phase diagram for 
model l|T|l: however, they have not investigated the uni- 
versality classes of the transitions in detail. We obtain 
the phase diagram here for density p = 1 by using the 
FSDMRG method which, as we show below, gives very 
accurate results for the nature of ordering in the differ- 
ent phases and the types and universality classes of the 
transitions. We restrict ourselves to the case of integer 
density (p — 1) since we want to explore the competition 
between SF, MDW, and MI phases. We note in passing 
that, even for V = 0, the Bose-Hubbard model (JTJ can- 
not be solved exactly unlike its fermionic counterpart; 
and for the fermionic case too there has been renewed 
interest in the phase diagram of the extended Hubbard 
model 2 1 . 

Before proceeding further we give a brief summary of 
our results. Our FSDMRG phase diagram for model QJ, 
with d = 1 and p = 1, is given in the (U,V) plane of 
Fig. (JTJ. It consists of three phases; Superfluid (SF), 
Mott Insulator (MI) and Mass Density Wave (MDW). 
For small values of the interactions U and V, the SF 
phase dominates as is to be expected since the bosons in- 
teract weakly here. However, as the interaction strengths 
increase, either MI or MDW phases get stabilized. The 
former dominates when U is much larger than V, whereas 
the latter dominates if U and V are both large and com- 
parable: A large, repulsive V disfavors a phase with a 
uniform density of bosons on nearest-neighbor sites; in- 
stead, an MDW phase, with a periodic variation of the 
boson density, is stabilized. The lattice we consider is 
bipartite and has two sublattices A and B (say odd- 
numbered and even- numbered sites); the ground state in 
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FIG. 1: The FSDMRG phase diagram of the one-dimensional, 
extended Bose-Hubbard model for density p — 1 showing SF, 
MI, and MDW phases. The two full lines indicate continuous 
transitions whereas the dashed line is a first-order bound- 
ary; these meet at a bicritical point. We do not consider the 
shaded region V > U. 

the MDW phase is, therefore, doubly degenerate since 
the peaks in the mass-density wave can lie either on 
the A or the B sublattice. If the bosons are charged 
this MDW phase is a charge-density-wave (CDW) phase. 
By using the FSDMRG method we have determined 
the phase boundaries between these phases. The MI- 
SF phase boundary in Fig. (JJJ lies in the Kosterlitz- 
Thouless (KT) universality class, whereas the MDW-SF 
phase boundary has both KT and two-dimensional-Ising 
characters. For large values of U and V, the MI-MDW 
transition occurs directly and is of first-order [dashed line 
in Fig. (I)]; as noted above, at t = 0, a first-order MI- 
MDW transition is obtained at U = 2V. Within the 
accuracy of our calculation, the MI-SF, MDW-SF, and 
MI-MDW phase boundaries meet at a bicritical point at 
(U, V) = (8.5±0.05,4.75±0.05). We have looked for, but 
not found, a supersolid phase with both SF and MDW 
order. A very brief discussion of some of our preliminary 
results has been given in Ref. pSjj. 

The remaining part of this paper is organized as fol- 
lows. Section II contains the details of our Finite-size 
Density-Matrix Renormalization Group (FSDMRG) cal- 
culation. Section III contains our results. We end with 
concluding remarks in Section IV. 

II. FSDMRG CALCULATIONS 

The Finite-Size Density-Matrix Renormalization 
Group (FSDMRG) method has proven to be very useful 



3 



in studies of one-dimensional quantum systemsL2Si2&. To 
make this paper self-contained we summarize the salient 
points of this method. Open boundary conditions are 
preferred for such calculations since the loss of accuracy 
with increasing system size is much less than in the 
case of periodic boundary conditions. The conventional 
FSDMRG method consists of the following two steps: 

1. The infinite-system density- matrix renormalization 
group method (DMRG), where we start with a sys- 
tem with four sites, add two sites at each step of 
the iteration, and continue till we obtain a system 
with the desired number L of sites ( in most of our 
calculations we use L ~ 100 but, in some represen- 
tative cases, we have gone up to L = 200). 

2. The finite-system method in which the system size 
L is held fixed, but the energy of a target state 
is improved iteratively by a sweeping procedure, 
described below, till convergence is obtained. 

For a model like Eq. we first construct the Hamilto- 
nian matrix of the superblock configuration • • BJ, 
where Bf and B^ represent left- and right-block Hamil- 
tonians, respectively, and each one of the • represents a 
single-site Hamiltonian. In the first step of the DMRG 
iteration both Bf and B^ also represent single sites, so, 
at this step, we have a four-site chain. We now diagonal- 
ize the Hamiltonian matrix of the superblock and obtain 
the energy and the eigenfunction of a target state. In 
our study the target state is the ground state of the sys- 
tem of size L with either N = L or N = L ± 1 bosons. 
The latter is required for obtaining the gap in the energy 
spectrum. We now divide the superblock into two equal 
halves, the left and the right parts, which are treated, re- 
spectively, as the system and the universe. The density 
matrix for this system, namely, B2 = B^ •, is calcu- 
lated from the target state. If we write the target state 
as I i>) = Hij^i,} I i) I i)> where | i) and | j) arc, 
respectively, the basis states of the system and the uni- 
verse, then the density matrix for the system has ele- 
ments pi^i = 'Ylij' l Pi,j^ l i',j- The eigenvalues of this den- 
sity matrix measure the weight of each of its eigenstates 
in the target state. The optimal states for describing 
the system are the ones with the largest eigenvalues of 
the associated density matrix. In this first step of the 
DMRG the superblock, and hence the dimension of the 
density matrix, is small, so all the states can be retained. 
However, in subsequent steps, when the sizes of the su- 
perblocks and density matrices increase, only the most 
significant states are retained, say the ones correspond- 
ing to the largest M eigenvalues of the density matrix (in 
our studies we choose M — 128). We then obtain the ef- 
fective Hamiltonian for the system Bj in the basis of the 
significant eigenstates of the density matrix; this is used 
in turn as the left block for the next DMRG iteration. 
In the same manner we obtain the effective Hamiltonian 
for the right block, i.e., B£ = »B^. In the next step 
of the DMRG we construct the Hamiltonian matrix for 



the superblock Bj • • Bj, so the size of the system in- 
creases from L = 4 to L = 6. For a system of size L, we 
continue, as in the first step, by diagonalizing the Hamil- 
tonian matrix for the configuration • •H r L _ and 

2 2 

setting Bi = B^_ • and B ^ = •H r L in the next step 
of the DMRG iteration. Thus at each step of the DMRG 
iteration the left and right blocks increase in length by 
one site and the total length L of the chain increases by 
2. 

In the infinite-system DMRG method outlined above 
the left- and right-block bases are not optimized in the 
following sense: The DMRG estimate for the target-state 
energy, at the step when the length of the system is L, 
is not as close to the exact value of the target-state en- 
ergy for this system size as it can be. It has been found 
that the FSDMRG method overcomes this problem^. In 
this method we first use the infinite-system DMRG iter- 
ations to build up the system to a certain desired size 
L. The L-site superblock configuration is now given by 

Bl . • • B r i _ 1 . In the next step of the FSDMRG 

2 1 2 1 

method, the superblock configuration B^ • • B^, , 

2 2 ~~ ^ 

which clearly keeps the system size fixed at L, is used. 
This step is called sweeping in the right direction since 
it increases (decreases) the size of the left (right) block 

by one site. For this superblock the system is B^», the 

2 

universe is •H r L , the associated density matrix can be 

found, and from its most significant states the new effec- 
tive Hamiltonian for the left block, with (-j + 1) sites, is 
obtained. We sweep again, in this way, to obtain a left 
block with (■§■ + 2) sites and so on till the left block has 
(L — 3) sites and the right block has 1 site so that, along 
with the two sites in between these blocks, the system 
still has size L; or, if a preassigned convergence crite- 
rion for the target-state energy is satisfied, this sweeping 
can be terminated earlier. Note that, in these sweeping 
steps, for the right block we need BJ to B£_ 3 , which 
we have already obtained in earlier steps of the infinite- 
system DMRG. Next we sweep leftward: the size of the 
left (right) block decreases (increases) by one site at each 
step. Furthermore, in each of the right- and left- sweep- 
ing steps, the energy of the target state decreases system- 
atically till it converges (we use a six-figure convergence 
criterion in our calculations). 

We use a slightly modified form of the FSDMRG 
method in which we sweep, as described above, at ev- 
ery step of the DMRG scheme and not only in the one 
that corresponds to the largest value of L. This helps 
in obtaining accurate /3 functions which we use to ob- 
tain critical exponents (see below) at continuous transi- 
tions. Furthermore, since the supcrfluid phase in mod- 
els such as Eq. Q, in d = 1 and at T — 0, is critical 
and has a correlation length that diverges with the sys- 
tem size L, finite-size effects must be removed by using 
finite-size scaling as we show below. For this purpose, 
the energies and correlation functions, obtained from a 
DMRG calculation, should have converged properly for 
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FIG. 2: The percentage difference between the ground state 
energies AE L = (eE mrg ~ E [ SDMRG )/E[ SDMRG obtained 
by DMRG and FSDMRG methods, for a representative point 
in the SF phase of Fig Q , plotted as a function of L. 



each system size L. It is important, therefore, that we use 
the FSDMRG method as opposed to the infinite-system 
DMRG method especially in the vicinities of continu- 
ous phase transitions. We find that convergence, to a 
specified accuracy for the target-state energy, is faster in 
the MI phase than in the SF phase. Figure shows, 
at a representative point in the SF phase of Fig. @J, 
how the percentage disagreement between DMRG and 
FSDMRG ground-state energies increases with L in our 
calculations. 

Since the bases of left- and right-block Hamiltonians 
are truncated by neglecting the eigenstates of the density 
matrix corresponding to small eigenvalues, this leads to 
truncation errors. If we retain M states, the density- 
matrix weight of the discarded states is Pm — — 
u> a ), where uj a are the eigenvalues of density matrix. Pm 
provides a convenient measure of the truncation errors. 
We find that these errors depend on the order-parameter 
correlation length in a phase. For a fixed M, we find 
very small truncation errors in the MI and MDW phases; 
these grow as the MI-SF and MDW-SF transitions are 
approached; and the truncation errors are largest in the 
SF phase. In our calculations we choose M such that the 
truncation error is always less than 5 x 10 -6 ; we find that 
M = 128 suffices. 

The number of possible states per site in the Bose- 
Hubbard model is infinite since there can be any number 
of bosons on a site. In a practical DMRG calculation 
we must restrict the number n max of states or bosons 
allowed per site. The smaller the interaction parameters 
U and V, the larger must n mnx be. As in earlier calcu- 



sufficient for the values of U and V considered here; we 
have checked in representative cases that our results do 
not change significantly if n max = 5. 

In summary, then, our FSDMRG procedure gives us 
the energy El (N) for the ground state of model (JTJ and 
the associated eigenstate | tPqln)- Given these we can 
calculate the energy gaps, order parameters, and corre- 
lation functions that characterize all the phases of this 
model and thence the phase diagram. We discuss this in 
Section III. 



III. RESULTS 

The single-particle energy gap Gl for a system of size 
L, the order parameter for the MDW phase, and the cor- 
relation functions that characterize SF and MDW phases 
in model Q can be defined in a straightforward manner 
in terms of the energies and wavefunctions mentioned in 
the previous Section. The energy gap is 

Gl=El{N + 1)+El{N~1)-2E l (N), (3) 

where El(N) is the ground-state energy for a system of 
size L with N bosons; since we are interested in studying 
the case p = 1, we increase the number of bosons by 2 
at every DMRG step in which 2 sites are added to the 
system (Sec. II) so that p = N/L = 1. We expect, and 
show explicitly below, that this gap is positive in both MI 
and MDW phases, which are incompressible insulators, 
but it vanishes in the SF phase, which is compressible. 

The correlation function that characterizes the SF 
phase is 



(r) = (ipoLN\ala r \ip LN), 



(4) 



where \iPoln) is the ground-state wavefunction of the sys- 
tem with size L and N bosons. The associated correlation 
length can be obtained from the second moment of this 
correlation function, namely, 



E,- 2 rf F (r) 



1/2 



(5) 



Note that ££ F is the correlation length for SF ordering 
in a system of size L; it remains finite so long as L < oo. 

The MDW phase can be differentiated from the MI 
phase by using the order parameter for the MDW phase 
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MDW 



is 



-l) l (-0OLJv|(n-i - p)\iPoln), (6) 



and the associated correlation function 



r MDW 
1 L 



(r) = (ipoLN\(n - p)(h r - p)\iPoln); (7) 



lation a 7 i 9 i 28 on related models, we find that 
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the correlation length for MDW ordering ^ L dDW can be 
defined as in Eq. © but with T A L IDW instead of Tf F . 

It is useful first to discuss the case V = 0. Here we re- 
produce the well-understood SF-MI transition. In d = 1 
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FIG. 3: L/£l F plotted as a function of U for different system 
sizes and V — in model Q. The coalescence of different 
curve for U ~ 3.4 shows a Kosterlitz-Thouless-type SF-MI 
transition. 



the appearance of the SF phase is signaled by the diver- 
gence of the correlation length £ l f j F ^ OCl . For a finite sys- 
tem ££ f is finite and we must extrapolate to the L — > oo 
limit, which is best done by using finite-size scaling. In 
the critical region the correlation length 



(8) 



where the scaling function f(x) ~ a; — > 0. Thus plots 
of L/t^if versus E7, for different system sizes L, consist 
of curves that intersect at the critical point at which the 
correlation length for L — oo diverges. Such plots are 
given in Fig. © for V = 0. Curves for different values of 
L coalesce for U < U c ~ 3.4 indicating the existence of 
a critical SF phase, with a diverging correlation length, 
i.e., a power-law decay of correlations, for all < U < U c . 
The single-particle gap Gl scales as the inverse of this 
correlation length and is, therefore, zero for L = oo for 
< U < U c . For U > U c we have an MI phase with a 
finite correlation length and a nonzero gap. Figure © 
also suggests that the SF-MI transition here is of the KT 
type. We can quantify this by calculating the /? function 
at this transition via the Roomany-Wyld (RW) approxi- 

mant9 7 i 29 4 30 
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He L F /( s L n/HL/L') 



(9) 



where L and L' denote two system sizes and £' L = 
d£,L F /dU. For a KT transition the correlation length 
Ckt ~ exp(c/(?7 - U c ) a ) and /3 ~ (U - U c ) 1+a , with 
a = 1/2. In Fig. Q we show the j3 function for the 
SF-MI transition for V = 0. To obtain this j3 func- 
tion we use L = 98 and L' = 100 in Eq. ©. Our 



FIG. 4: The f3 function for the MI-SF transition for V = 
obtained by using RW approximants [Eq. JSJ] with L — 98 
and L' = 100. The full line is a fit to the form 0(U) = 
c(U - t/ c ) (1+CT) ; we get c = 0.37 ± 0.01, U c = 3.35 ± 0.02 and 
a = 0.48 ±0.05. 



fit to the data of Fig. gj yields U c = 3.35 ± 0.02 and 
a = 0.48 ± 0.05 which are consistent with the values re- 
ported earlier 7-9 . If we fit our data (here and below) over 
a fixed region of (U — U c )/U c , then our nonlinear least- 
squares programme yields smaller errors. The conserva- 
tive error bars we quote here reflect the range over which 
our fitted values for U c , cr, etc., vary when we change the 
region of (U — U c )/U c over which we fit our data, namely, 
0.001 < \U-U C \/U C < 0.2 or 0.001 < \U-U C \/U C < 0.35. 

Initially the nearest-neighbor interaction V suppresses 
the SF phase relative to the MI phase, but at larger val- 
ues of V this trend is reversed leading to a reentrant SF 
[Fig. 0]. Figure (JSJi) shows a plot of L/^ F versus U 
for V — 0.5; by comparing this with Fig. we see that 
U C (V = 0.5) < U C (V = 0). Curves for different values 
of L coalesce for U < U C (V — 0.5) ~ 3.0 indicating a 
power-law SF phase, as for V = 0, and a KT-type MI-SF 
transition. Again we use RW approximants to obtain the 
(3 function shown in Fig. © and our numerical fit yields 
U c = 2.95 ± 0.02 and cr = 0.53 ± 0.05. Eventually the 
MI-SF phase boundary in Fig. Q turns back; a represen- 
tative plot oiL/^ff versus L, for V = 2.5, illustrates this 
[Fig. J3>)]; we find U C (V = 2.5) ~ 4.2. This reentrance 
of the SF phase, with increasing V, was not resolved by 
the study of Ref. 0. 

For sufficiently large values of V we can have an MDW 
phase and an SF-MDW transition, at small values of U, 
and an MI-MDW transition at large values of U as shown 
in the phase diagram of Fig. QJ. In Fig. Q we plot the 
MDW order parameter Mmdw as a function of 1/L for 
U = 6 and values of V ranging from V = 3.0 to 5.3 
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FIG. 5: L/f£ F plotted as a function of U for different system 
sizes and for (a) V = 0.5 and (b) V = 2.5 in model JIJ. The 
coalescence of different curve for U < 3.0 in (a) and U < 4.2 
in (b) show KT-type SF-MI transitions. 
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FIG. 6: The f3 function of the MI-SF transition for V = 0.5, 
obtained by using the RW approximants with L — 98 and 
L' = 100 as in Fig. @. The full line is a fit to the form 
(3{U) = c(C/-C/ c ) (1+CT) ; we get c = 0.27±0.01, U c = 2.95±0.02 
and a = 0.53 ± 0.05. 



in steps of 0.1. We see that Mmdw goes to zero for 
V < V C (U = 6) ~ 3.9 whereas it is nonzero for higher 
values of V. To determine the universality classes of 
the MI-SF and SF-MDW transitions we plot in Fig. H 
both the order parameter Mmdw an d L/^ff as func- 
tions of V for U = 6. The different curves for L/^f F 



FIG. 7: The MDW order parameter Mmdw versus 1/L for 
U — 6 and different values of V starting from 3.0 below to 5.3 
above in steps of 0.1. 



coalesce in the region 3.5 < U < 3.9 indicating an SF 
phase sandwiched between MI and MDW phases. Both 
SF-MI and SF-MDW transitions are continuous. To con- 
firm this we have also obtained plots of the ground-state 
energy E = Huil-^oo El{N) as a function of V for fixed 
U. In Fig. @ we plot E and dE /dV versus V; this 
plot shows no discontinuity at the SF-MI and SF-MDW 
transitions (as it does at the first-order MI-MDW tran- 
sition discussed below for U = 12). To determine the 
universality classes of these SF-MI and SF-MDW tran- 
sitions we have obtained (3 functions, via RW approxi- 
mants for £f F , in Fig. for C7 = 6.0. For the MI-SF 
transition, we get V c = 3.59 ± 0.05 and a = 0.47 ± 0.05; 
and, for the SF-MDW transition, V c = 3.78 ± 0.05 and 
(7 = 0.49 ± 0.05. Thus both of these transitions are of 
the KT type. However, in addition, the MDW-SF tran- 
sition also has an Ising character (two-dimensional) since 
the MDW phase has a doubly degenerate ground state 
as mentioned above. To extract such Ising-model ex- 
ponents, we use the form Mmdw ~ (V — V C Y MDW as 
V I V c , where Pmdw is the MDW order-parameter ex- 
ponent. For U = 6 our fit yield V c = 3.87 ± 0.05 and 
Pmdw = 0.12±0.01 (Fig.HT)l in good agreement with the 
two-dimensional, Ising order parameter exponent. Note 
that the value of V c obtained from this fit for Mmdw is 
within error bars of that obtained from the [3 function for 
the SF-MDW transition. Thus, within our calculation we 
cannot resolve a supersolid phase which has long-range 
SF correlations and MDW ordering. 

From the phase diagram of Fig. (JIJ we see that, for suf- 
ficiently large values of U, there is no SF phase and only a 
direct, first-order MI-MDW transition. This direct tran- 
sition shows up clearly in Fig. 1121 where we plot L/^f F 
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FIG. 8: Plots of L/Cf F and Mmdw versus V for U = 6 and 
L — 90 and 100. The coalescence of different curves of F , 
for 3.5 < V < 3.9, shows an SF phase sandwiched between 
MI and MDW phases; for 3.9 < V we obtain the MDW phase 
with Mmdw > 0. 
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FIG. 9: Plots of the ground-state energy Eq and its derivative 
dEo/dV versus V for [7 = 6. There is no visible discontinu- 
ity in dEo/dV so both MI-SF and SF-MDW transitions are 
continuous. 



and Mmdw versus V for U = 12: Curves of L/^jf, 
for different values of L, do not merge at any point, so 
we can conclude that no power-law SF phase intervenes 
between MI and MDW phases. Furthermore, the sharp 
jump in Mmdw at V ~ 6.3 indicates that we have a 
first-order MI-MDW transition. This is corroborated by 
the plots of the ground-state energy E Q and its derivative 
dE/dV given in Fig.^|for U = 12; the discontinuity of 




FIG. 10: The j3 functions for MI-SF (a) and MDW-SF (b) 
transition for [7 = 6, obtained by using RW approximants 
with L = 98 and U = 100. The full line is a fit to the form 
P(U) = c{U - Uc) (1+rT) . For the MI-SF transition (a) we 
get c = 0.76 ± 0.01, U c = 3.59 ± 0.05 and a = 0.47 ± 0.05 
and for the MDW-SF transition (b) we get c = 0.86 ± 0.01, 
U c = 3.78 ± 0.05 and a = 0.49 ± 0.05. 
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FIG. 11: A log-log plot (base 10) of the MDW order parame- 
ter Mmdw versus V — V c for U = 6. The straight line is a fit 
to the form Mmdw ~ (V - V C ) 0MDW , V I V c . Our fit yields 
Vc = 3.87 ± 0.05 and [5m dw = 0.12 ± 0.01. 
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FIG. 12: Plots of LUl F and Mmdw versus V for U = 12 and 
different system sizes L. Note that the curves for L/£l do not 
meet at any point; however, Mmdw jumps at the first-order, 
MI-MDW transition at V ~ 6.3. 
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FIG. 14: Plots of L/Cf versus V for different values of U and 
different system sizes L(80, 90, 100). The region over which 
curves for different values of L coalesce decreases as we go 
from U = 8 to C/ = 8.5; and there is no coalescence for U = 9. 
This shows how the SF phase shrinks as we approach the 
bicritical point. 



region over which the curves of L/^f j F coalesce) shrinks 
as we approach the bicritical point. Similarly Fig. Q15(l 
shows how the jump in Mmdw across the MI-MDW first- 
order transition decreases as we approach the bicritical 
point. To confirm that this is indeed the topology of 
the phase diagram in the region where the MI-SF, SF- 
MDW, and MI-MDW phase boundaries meet, we must 
obtain the critical exponents in the vicinity of the bicriti- 
cal point. This is beyond the accuracy of our calculation 
at the moment. Nor can we rule out completely more 
complicated topologies of phase diagrams in which very 
closely spaced tricrital points and critical endpoints are 
used to link the three phase boundaries we have studied 
above. 



FIG. 13: The ground-state energy Eo and its derivative 
dEo/dV versus V for U = 12.0; the jump in dEo/dV at 
V ~ 6.3 shows that the MI-MDW transition is first order. 

dE/dV at V ~ 6.3 indicates the first-order nature of the 
transition. 

The only feature of the phase diagram of Fig. QJ that 
remains unexplored now is the region in which the con- 
tinuous MI-SF and SF-MDW phase boundaries meet the 
first-order MI-MDW boundary. The simplest topology 
possible here is that these meet at one bicritical point. 
Our data are not inconsistent with such a topology. Fig- 
ure (|14f> shows, via plots of L/££ F versus V, for different 
values of L and U, how the extent of the SF phase (the 



IV. CONCLUSIONS 

We have carried out an extensive study of the one- 
dimensional, extended Bose-Hubbard model QJ by using 
the FSDMRG method. Our study yields ground-state 
energies, single-particle gaps, the MDW order parame- 
ter, and SF correlation functions and correlation lengths. 
By studying these we obtain an accurate phase diagram 
[Fig. 0J] for this model. This shows continuous MI- 
SF and SF-MDW transitions meeting the first-order MI- 
MDW boundary at a bicritical point. 

We have looked for, but not found, a supersolid (SS) 
phase which, in the context of the lattice model we study 
here, would exhibit power-law superfluid correlations, as 
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FIG. 15: The MDW order parameter Mmdw plotted ver- 
sus V/V c for different values of U across the MI-MDW first- 
order transition. The jump in Mmdw at the first-order phase 
boundary decreases as we approach the bicritical point at 
(U, V) = (8.5 ± 0.05, 4.75 ± 0.05). 



in the SF phase, and a nonzero order parameter Mmdw, 
as in the MDW phase. It is likely^i that further than 
nearest-neighbor interactions will be required to stabilize 
the SS phase as we will explore elsewhere. 

We hope our detailed study of model (Q) will stimulate 
experimental studies. A recent experimental studySS has 
shown that it is possible, by a suitable choice of con- 
fining potentials in optical lattices, to obtain a physical 
realization of the one-dimensional Bose-Hubbard models. 
It would be interesting to see how nearest-neighbor inter- 
actions, like V, can be obtained in such lattices. If this 
can be done, the rich phase diagram of Fig. can be 
explored experimentally. 
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